%% 通用变量
ku=10; % 流入阀
 
%% 单容水槽
k=2;% 这个是流阻相关
H0=5;% 水槽初始高度
F=20; %面积
R=2*sqrt(H0)/k;
K=ku*R;
T=R*F;

%% 双容水槽
H10=5; %水槽1初始
H20=1; %水槽2初始
k1=2;
k2=3;
F1=20;
F2=10;
R1=2*sqrt(H10)/k;
R2=2*sqrt(H20)/k;
T1=R1*F1;
T2=R2*F2;
K1=ku*R1;
K2=R2/R1;
r=R2/(R1+R2);
K_shuang=r*ku*R2;

% H20_1=5
H20_1=5;
R1_1=2*sqrt(H10)/k;
R2_1=2*sqrt(H20_1)/k;
T1_1=R1_1*F1;
T2_1=R2_1*F2;
K1_1=ku*R1_1;
K2_1=R2_1/R1_1;
r_1=R2_1/(R1_1+R2_1);
K_shuang_1=r_1*ku*R2_1;

% H20_2=10
H20_2=10;
R1_2=2*sqrt(H10)/k;
R2_2=2*sqrt(H20_2)/k;
T1_2=R1_2*F1;
T2_2=R2_2*F2;
K1_2=ku*R1_2;
K2_2=R2_2/R1_2;
r_2=R2_2/(R1_2+R2_2);
K_shuang_2=r_2*ku*R2_2;

% H20_3=15
H20_3=15;
R1_3=2*sqrt(H10)/k;
R2_3=2*sqrt(H20_3)/k;
T1_3=R1_3*F1;
T2_3=R2_3*F2;
K1_3=ku*R1_3;
K2_3=R2_3/R1_3;
r_3=R2_3/(R1_3+R2_3);
K_shuang_3=r_3*ku*R2_3;

% H20_4=20
H20_4=20;
R1_4=2*sqrt(H10)/k;
R2_4=2*sqrt(H20_4)/k;
T1_4=R1_4*F1;
T2_4=R2_4*F2;
K1_4=ku*R1_4;
K2_4=R2_4/R1_4;
r_4=R2_4/(R1_4+R2_4);
K_shuang_4=r_3*ku*R2_4;